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Recent theory and experiment have revealed that strong spin-orbit coupling can have dramatic 
qualitative effects on the band structure of weakly interacting solids. Indeed, it leads to a distinct 
phase of matter, the topological band insulator. In this paper, we consider the combined effects of 
spin-orbit coupling and strong electron correlation, and show that the former has both quantitative 
and qualitative effects upon the correlation-driven Mott transition. As a specific example we take Ir- 
based pyrochlores, where the subsystem of Ir 5d electrons is known to undergo a Mott transition. At 
weak electron-electron interaction, we predict that Ir electrons are in a metallic phase at weak spin- 
orbit interaction, and in a topological band insulator phase at strong spin-orbit interaction. Very 
generally, we show that with increasing strength of the electron-electron interaction, the effective 
spin-orbit coupling is enhanced, increasing the domain of the topological band insulator. Further- 
more, in our model, we argue that with increasing interactions, the topological band insulator is 
transformed into a "topological Mott insulator" phase, which is characterized by gapless surface 
spin-only excitations. The full phase diagram also includes a narrow region of gapless Mott insu- 
lator with a spinon Fermi surface, and a magnetically ordered state at still larger electron-electron 
interaction. 

PACS numbers: 71.30.+h, 71.70.Ej, 73.20.At 

Introduction: The spin-orbit interaction (SOI), though apparently a "weak" relativistic correction to the Schrodingcr 
equation (outside of high energy physics), is coming increasingly to the fore in modern condensed matter physics. The 
discovery of Topological Band Insulators (TBIs) in theory [1, 2, 3, 4] and experiment [5, 6] has revealed a surprising 
omission in the "textbook" Bloch theory of the electronic structure of weakly correlated solids. In these remarkable 
materials, strong spin-orbit interactions allow a non-trivial topology of the electron bands, resulting in protected 
"helical" edge and surface states in two and three dimensional systems. Many other interesting phenomena, including 
quantum number fractionalization and magneto-electric effects have been predicted to occur in these systems, and are 
the subjects of a growing experimental effort. In parallel, strong SOIs have been identified in a growing variety of Mott 
insulators, in which the insulating behavior is driven by electron correlation rather than band structure. For instance, 
SOIs arc likely responsible for the large Wilson ratios observed in many frustrated magnets at low temperature [7, 8], 
and may be the driving force for the formation of a "spin-orbital liquid" in some Fe-spinels [9]. They have been 
experimentally shown to control the orbital state in the Ir oxide Sr 2 Ir04 using resonant x-ray scattering [10]. A 
natural question is how these two classes of phenomena are connected - how does a material progress from weak to 
strong correlation with strong SOIs? This is the subject of the Mott transition with strong SOIs. 

In the search for strong SOIs, one is driven to consider materials based on heavy atoms (as SOI scales with atomic 
number as Z 4 ). So far, the experimental examples of topological insulators have been based on atoms (Hg and 
Bi) in which the active electrons closest to the Fermi energy inhabit s and p orbitals, which are generally weakly 
correlated. However, strong spin-orbit coupling and electron correlation may be expected in compounds involving d 
electrons. Especially in 5d transition metals, correlation effects are often of weak or intermediate strength, making 
them perhaps competitive with those of SOIs and band structure. Due to the surge of activity in Ir-based transition 
metal oxides [7, 10, 11, 12, 13, 14, 15], we choose one of these as our model system, and study a simple phenomenological 
Hamiltonian for the Mott transition in Ir pyrochlores. We will see that this crystal structure naturally fosters TBI 
physics in the weak-correlation limit, and therefore allows for a study of the most interesting regimes of the Mott 
transition with strong SOIs. 

Our main results are as follows. Most importantly, we find that Mott correlations enhance the effects of SOIs, 
due to the suppression of effective bandwidth. This leads to increased stability of the topological insulator state 
with interactions. In addition, we suggest the existence of an intermediate Topological Mott Insulator (TMI) state, 
between the TBI and the familiar Mott insulating phases at stronger interactions. The TMI phase may be regarded 
as a quantum spin liquid state, in which the electrons are effectively spin-charge separated, and neutral fermionic 
spin-1/2 spinons are deconfincd. The non-trivial band topology of the TBI is transferred to the spinons, which are 
thereby gapped in the bulk but gapless at the surface. These results are summarized in the phase diagram in Fig. 1. 

Model: We begin by constructing a minimal microscopic model for the Mott transition in Ir-based pyrochlore oxides, 
A 2 Ir 2 07, where A is assumed to be an inert ion in the 3+ state. Therefore we consider Ir 4+ ions on the pyrochlore 
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FIG. 1: Phase diagram based on the slave rotor approximation and strong coupling limit, as a function of Hubbard repulsion 
U and spin-orbit coupling A (relative to hopping t). The four main phases occuring for moderately strong electron-electron 
repulsion are a Metallic phase, Topological Band Insulator (TBI), Topological Mott Insulator phase (TMI), and Gapless Mott 
Insulator (GMI). The dashed line denotes an additional zero-gap semiconductor state due to an "accidental" gap closing. The 
dotted line schematically separates the large-[/ region, where magnetic ordering is expected. As discussed in the main text, 
long-range Coulomb interactions are expected to induce an excitonic region in the vicinity of the Mctal-TBI boundary shown 
here. 



lattice, which is viewed as the FCC lattice with a tetrahedral basis, Fig. 2(a). Each Ir ion is surrounded by an oxygen 
octahedron, Fig. 2(b). In general, these octahedra are somewhat distorted, dependent upon the difference of the 
oxygen displacement parameter x from the ideal value of 5/16. We neglect this small distortion (of 5 and 8 percent 
for A=Pr,Eu) here, in which case the Ir electrons experience approximately cubic Oh symmetry, and the crystal field 
splits the five d-orbital states into (upper) e g doublet and (lower) t 2g triplet. The value of the cubic crystal field 
splitting is assumed large compared to the Hund exchange energy, thus the latter is neglected. The on-site spin-orbit 
interaction thus is taken to act within the t 2g manifold. In local axes bound to the oxygen octahedron the tig orbitals 
are d yz ,d zx and d xy . In Supplementary Information we present the rotation matrices that turn the cubic axes into 
the local octahedral ones at each Ir site. The axes themselves are not shown in Fig. 2 for clarity reasons. 



FIG. 2: a)Pyrochlore lattice of corner-sharing tetrahedra. It can be viewed as the FCC lattice with tetrahedral bases added at 
each site. One of such bases, with Ir ions numbered from to 3, is shown, b) Oxygen mediated hopping between Ir sites. Sites 
(on the left) and 3 (on the right) of the tetrahedral basis are shown (large gray spheres), together with their oxygen octahedral 
environment (small red spheres). On the "shared" oxygen site we show its p y orbital (green) with respect to the coordinate 
system of site 0, and p x orbital (blue) with respect to the coordinate system of site 3. Belonging to different coordinate systems, 
these orbitals are not orthogonal, the angle between them being « 84°. Electrons can hop from local d yz orbital on site onto 
p y , and from the d xy orbital on site 3 onto p x . Since the two p orbitals are not orthogonal, an effective Ir-Ir hopping is induced. 

Now consider the SOI, denoting its strength by A. It is well-known that, projected within the tig manifold, the 
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FIG. 3: Electronic band structure of Ir 5d electrons on the pyrochlore lattice at large spin-orbit coupling, X/t — > oo. Only the 
relevant four doubly degenerate bands are shown. A band gap between the filled lower two bands and empty upper two bands 
is clearly seen. 



orbital angular momentum operator is equal to minus the usual spin-one angular momentum operator L [16]. Thus we 
may effectively consider the pseudo-total angular momentum J — L+ S, which commutes with the SOI. The multiplet 
structure for the usual £ = 1 states is inverted, and the t 2g manifold splits into a j = 1/2 doublet with energy A, and 
a j = 3/2 quadruplet with energy —A/2. We will work in the local diagonal basis of the j eigenstates, and introduce a 
single label a such that a = 1, 2 and a = 3 ... 6 denote the doublet and quadruplet, respectively, with orbital energy 
e a = A for the doublet and e a = —A/2 for the quadruplet. The specific form of the wave functions of these states are 
given in the Supplementary Information. 

We assume that hopping between nearest neighbor Ir ions is accomplished via the oxygen atoms nearest to a given 
pair, Fig. 2(b). In reality, this is not necessarily the case [17]. However, our assumption minimizes the number of 
free parameters, and is resilient to perturbations that are not too strong. The model with oxygen-mediated hopping 
is preferable as it contains a single parameter determining the hopping strength: the hopping integral between Ir t 2g 
states and O p-orbitals (V p d^ in the terminology of Ref. [18]). Integrating out the oxygens, and taking the simplest 
on-site Coulomb interaction involving the total charge only, we arrive at the Hubbard Hamiltonian, 

H = E^-^L^ffia + i E T aL'd R ^d m>a , + jJ2[J2 d ma d ma-^d) , (1) 

Ria [R:i,R'i') Ri \ a J 

OLOL 

where R and i label the sites of the Bravais lattice and the tetrahedral basis, and rid = 5 is the number of 5d-electrons 
on each Ir 4+ ion. The strength of the hopping is parameterized by a single energy scale t oc V^ 2 d7r /A, where A is the 

energy difference between the Ir d and O p states. The procedure to obtain the dimensionless hopping matrices T™ a ,, 
arising from the Ir-OTr hopping path, taking into account the rotation between the local cubic axes of each Ir ion, is 
given in the Supplementary Information. 

Band structure: The Hamiltonian (1) contains two dimensionless parameters: X/t and U/t, which define the phase 
diagram in Fig. 1. It is instructive to consider first various simple limits. For U = 0, we have a free electron model, 
which is of course exactly soluble. Due to inversion symmetry, one obtains in general 12 doubly degenerate bands. 
For small X/t, these overlap at the Fermi energy and one obtains a metal. For large X/t, the upper 4 bands originating 
from the j — 1/2 doublet become well-separated from the lower 8 bands. Because there are four holes per unit cell, 
the upper 4 bands are half-filled in total. On inspection, we see (Fig. 3) that they exhibit a band gap, indicating the 
formation of a band insulator at large A. 

As shown by Fu and Kane [19], one can determine the band topology of an insulator with inversion symmetry cither 
from the parity of the Hamiltonian eigenstates at time-reversal invariant momenta, or from the number of Dirac points 
on the surface of the insulator. Applying the first criterion (see Supplementary Information), we find that the large 
X/t state is a pure "strong" TBI of the spinons (the weak Z 2 invariants vanish, consistent with cubic symmetry). We 
also calculated the surface state spectrum (Fig. 4), which shows the required odd number of intersections with the 
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FIG. 4: Surface state spectrum. Projected bulk spectrum (gray), and the surface spectrum (solid blue line) for a (010) surface 
are shown. The spectrum is plotted as a function of the surface wave vector k sur f = {k x ,k z ) running from T — (0,0) to 
X=(0,27r). 



Fermi level on passing between time reversal invariant surface momenta. The behavior for general X/t is as follows. 
For X/t < 2.8, one obtains a metallic state, while for X/t > 2.8 the bands separate at the Fermi energy. For almost 
all of this range of large X/t, the system is a (strong) TBI. However, an "accidental" closing of the band gap occurs 
at X/t f=a 3.3, at which point it is a zero-gap semiconductor with 8 Dirac points located along the (111) directions in 
reciprocal space. Because the number of these Dirac points is even, there is no change in band topology due to the 
gap closure. 

Strong coupling limit: Now consider large U/t. In this limit, one has a Mott insulator, and the Hamiltonian 
is effectively projected into the space of one hole per Ir site, and the system is described by a spin-orbit (Kugel- 
Khomskii type) model. Superexchange leads to spin (and orbital) exchange of order J ~ t 2 /U <C t. We see that 
in this limit SOI is only weak if A <C J <C t. Thus the strong SOI regime is greatly enhanced with increasing 
correlations, as the relevant "bandwidth" for large U is exchange rather than hopping. Since the general spin-orbital 
Hamiltonian is cumbersome, and only relevant for very weak SOI, we will focus only on the strong SOI regime. Here, 
only the half- filled doublet at each site is relevant, and the effective Hamiltonian is of Heisenberg spin-exchange type, 
with an effective spin 1/2 at each site. The derivation of the low-energy Hamiltonian by second order degenerate 
perturbation theory is standard. It is customary to write the resulting Hamiltonian as a sum of isotropic exchange, 
Dzyaloshinskii-Moriya (DM), and anisotropic exchange parts (suppressing the Bravais lattice index R for brevity): 
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The convention to make the choice of Da> unique is that i! > i. Since all -D's and T's are related by symmetry, it is 
enough to specify them for one bond only: 
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We note that the anisotropic exchange term, F qi, is small compared to the other ones. Spin ordering in this model 
was considered in Ref. [20]. Our DM term corresponds to the "indirect" case of Ref. [20], and \D i\/J as 0.63, i.e. 
very large DM interactions. In this case, a magnetically ordered ground state is expected, which breaks point group 
symmetries but does not enlarge the unit cell. 
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Slave-Rotor approach: An exact or accurate numerical solution for the full phase diagram for Eq. (1) is very 
challenging, especially at intermediate U/t. To study it, we employ the slave-rotor approximation of Ref. [21]. This 
approach has a number of merits. It becomes exact for U/t = 0, and captures the bandwidth reduction with increasing 
U/t. Its predictions for Hubbard models on other frustrated lattices without SOI at intermediate U /t are in agreement 
with more controlled approaches such as the path integral rcnormalization group [22], Gutzwiller-type variational 
wavefunctions [23], and the variational cluster method [24]. Since we will see that the Mott transition occurs at 
smaller U/t with increasing X/t, we expect that the slave-rotor approximation should be reasonable to describe it for 
the full range of SOI. It clearly fails at large U/t, but we can substitute direct analysis of the spin-orbital model in 
that limit. 

We decompose the physical electron annihilation operator as dm a — e~ l9fti fm a , where the angle 9m is the con- 
jugated variable to the number of electrons on site R, i (the "angular momentum" of the rotor), and the "spinon" 
fma carries the rest of the degrees of freedom. The constraint L Ri = J2 a fmafma — n di restricting the physical 
part of the Hilbcrt space, is treated on average. Further, we use the mean-field decomposition of the hopping term 
that couples the spinons and rotors according to AB — > A (B) + B (A) . This Mean Field Theory (MFT) reduces the 
Hamiltonian (1) to two uncoupled Hamiltonians for spinons and rotors: 

(m,R f i f ) 

aa 

He = ^^2L 2 Ri + h{L m + n d )+tQ e £ e **«-«W . (4) 
m (m,R'i') 

Here Lju = ~i-Q§^:, the coordinate-independent Lagrange multiplier h is introduced to treat the constraint on the 
angular momentum, and the couplings Qf and Qg need to be determined self-consistently from Qf = (e 19 "^ 10 ^') , 
Qe = J2 aa ' Ta a > ^/flj a /fli'a' ) ( n °tc Qe < in the self-consistent solution). Here we have made the so-called "uniform" 
mean-field approximation, on the grounds that it is the one which smoothly connects to the U/t = limit, and hence 
should be appropriate for small to intermediate U/t, where, anyway, we apply the approximation. Note that the 
strength of the spin-orbit interaction, A, is not renormalized. After the mean-field decomposition, the spinon and 
rotor sectors can be solved independently. 

We first consider the spinon Hamiltonian, Hf. It is identical to the non-interacting electron Hamiltonian, but with 
renormalized hopping t c g = tQf. All the preceding analysis carried out for U = can therefore be carried over with 
this replacement. Since Qf < 1, we indeed observe the dimcnsionlcss spin-orbit strength X/(Qft) is enhanced by 
correlations. Physically, however, we must take care that the / fermions are spinons and therefore their properties do 
not necessarily translate directly to the physical electrons. 

We now turn to the charge (rotor) sector. Hg describes rotor bosons moving on the pyrochlore lattice. The 
parameter h must be fixed by charge neutrality, {fjn a fm a ) = n <ii an d hence (Lin) — 0. We therefore take h = 0, 
which guarantees the latter condition, since then Hg has particle-hole symmetry, L Ri — > —LruOri — * —0Ri- Hg is 
then expected to exhibit two phases. For U/(Qgt) <C 1, the rotors are condensed, (e t6Ri ) ^ 0, while for U/(tQg) ^> 1, 
they form an uncondenscd Bose Mott insulator with a gap and (e Ri ) = 0. In the rotor condensate, the electron 
and spinon operators are proportional, and one has "Fermi liquid" states where the spinon band structure describes 
physical electrons. When the rotor is uncondensed, there is spin-charge separation and the spinons are emergent 
neutral quasiparticles carrying spin only. Separating these two states is expected to be a continuous Quantum 
Critical Point (QCP). 

Combining the solutions of the spin and charge sectors, one finds four possible phases. When the rotors are 
condensed, we obtain metallic and TBI phases. When the rotors are uncondensed, the corresponding phases are 
"exotic" non-magnetic insulators. For small X/(Qft), this is a gapless insulator with a spinon Fermi surface. For large 
X/(Qft) there is a bulk gap, but the spinons have the non-trivial TBI topolgy; we call this the "Topological Mott 
Insulator" (TMI) phase. Its properties will be discussed further below. 

To determine the boundaries between these phases, we must make MFT self-consistent. From the solution of Hf, 
we can obtain one self-consistent condition Qg = f{-^j), with /(£) determined numerically. The solution of Hg in 

turn fixes the other self-consistent equation, Qf — ff(^~t)- I n principle, <?(£) (which decreases with increasing £ from 
1 to 0) could be determined numerically, since Hg is amenable to quantum Monte Carlo without a sign problem. 
However, we content ourselves with a simple analytical approximation which produces reasonable results. Further 
details on the selfconsistent solution of MFT to extract the phase boundaries X(U) are given in the Supplementary 
Information. Here we present the resulting phase diagram only, Fig. 1. 
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Discussion: The combined predictions of the slave-rotor theory and the strong coupling limit are shown in the 
phase diagram of Fig. 1. Focusing first on weak to intermediate correlation, the theory predicts two types of metal- 
insulator transitions. The conventional Mott transition in which the insulator is of Mott type - the Gapless Mott 
Insulator (GMI) - occurs only for rather weak SOI. Instead, the dominant metal-insulator transition is actually to 
TBI rather than a Mott insulator, due to the bandwidth reduction with increasing U. This is not a Mott transition 
per se, as the TBI is qualitatively described by band theory. Instead, the true Mott transition occurs within the bulk 
insulating state, between the TBI and the Topological Mott Insulator (TMI). In the TMI, the charge is stripped from 
the spinons, resulting most notably in a change from a conducting surface to an insulating one, which nevertheless 
still supports gapless spin excitations. 

The above slave-rotor MFT, by construction, leads to exotic Mott insulators (the GMI, TMI) with spin-charge 
separation. More precisely, these states are U(l) spin liquids, which, when fluctuations are properly taken into 
account, must support an emergent dynamical U(l) gauge field. Such spin liquid states are known to be stable 
in three dimensions. In addition to the excitations already described, they support gapped "magnetic monopolcs", 
which are sources of the emergent magnetic flux. Such a monopolc in a topological insulator carries a quantized gauge 
charge [25, 26, 27, 28]. Deep in the TMI, where charge excitations are costly, we thus expect a spinon to be bound 
to the monopole, i.e. the monopole has spin. The low energy effects of gauge fluctuations in the GMI (and near 
the metal-GMI transition) have been discussed elsewhere (see Ref. [29]). We expect only moderate effects of gauge 
fluctuations on the gapless surface spinons in the TMI state, owing to the three-dimensional character of the gauge 
field. 

Is the mean-field starting point itself reasonable? Because in these phases the spinons behave very similarly to 
electrons in the nearby (metallic and TBI) conventional states, they are most natural in "weak" Mott insulators, in 
which there are significant charge fluctuations. This is in line with the MFT prediction of continuous Mott transitions 
for these phases. In reality, many effects beyond MFT might lead to these transitions becoming first-order, as often 
observed in experiment. However, numerous recent analytical and numerical studies [22, 30, 31, 32, 33, 34] support 
the existence of continuous or only weakly discontinuous Mott transitions in frustrated geometries such as the one 
studied here. Because the Mott transition is suppressed to small U /t by increasing SOI, we believe the predictions of 
the slave-rotor MFT are reasonable both for the GMI and the TMI phases. 

Turning specifically to the character of the weak to intermediate coupling transitions in the phase diagram, we 
expect that effects beyond MFT will indeed make several of them weakly first order. First, the metal-TBI transition, 
which in MFT is described by a simple crossing of valence and conduction bands at the V and L points in the Brillouin 
zone, respectively, is modified by the effects of long-range Coulomb interactions. This leads to exciton condensation, 
occupying a region of the phase diagram in the vicinity of the putative metal-TBI transition [35] . If electron-phonon 
coupling is not too strong, the most likely excitonic condensation is predominantly triplet, and this state should 
evince weak magnetic order [35]. Dependent upon the details of the magnetic order, band structure, and disorder, 
the transitions to this state might be first or second order. Second, the TBI-TMI transition, described in MFT by the 
condensation of the rotor, is modified by coupling to the emergent U(l) gauge field in the TMI state. This problem 
maps onto the abelian Higgs transition in 3+1=4 dimensions, where it can be studied by renormalization group 
methods. The gauge fluctuations destabilize the continuous mean-field transition, rendering it weakly first order [36]. 
The metal-GMI transition has been studied recently in Ref. [29], and is potentially continuous even with fluctuations. 
We do not discuss the GMI-TMI transition here. 

The slave rotor approximation neglects possible magnetic instabilities, which become increasingly prominent with 
increasing U/t. For comparison, in Ref. [17], a phase diagram in the same parameter space was obtained for the 
layered material Na 2 Ir03, based on density functional theory, in which the Mott phases are all magnetic. None of 
the phases we have found are perturbatively unstable to magnetic order (excepting the excitonic region described 
above), so the simplest scenario is one in which a direct transition occurs at larger U/t from the GMI/TMI states to 
magnetically ordered ones (e.g. along the dotted line in Fig. 1). It will be interesting to consider whether this might 
occur by condensation of the monopole excitations. It is also possible for magnetism to occur at smaller U /t for other 
reasons. For instance, nesting in the metallic state would promote spin density wave (SDW) formation, which could 
occur by a Slater transition below the Mott transition predicted by the slave rotor approximation. In our calculations, 
we do not see obviously nesting in the metallic Fermi surfaces, but this could occur if other hopping processes were 
important. Another interesting possibility is emergence of magnetic ordering at the surface of the TBI or TMI states, 
where the surface fermions can undergo Slater transitions. 

We conclude with a discussion of experiments and materials. Experiments on Ir pyrochlores, Ln2lr207 show a zero- 
temperature metal insulator transition with increasing Ln radius, occuring between Ln=Nd and Ln=Pr [12]. The 
insulating materials show a weak continuous T > metal-insulator transition, with some magnetic order evident in 
the low-temperature insulating state. This might be consistent with the excitonic state expected near the metal-TBI 
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boundary. However, numerous other interpretations are possible, and clearly more theoretical and experimental work 
is required to better understand these materials. More generally, the form of the phase diagram in Fig. 1 may apply to 
a variety of 4d and 5d transition metal compounds in which both spin-orbit coupling and correlations play important 
roles, and whose lattice structures are favorable for TBI physics in the weak correlation limit. In Ir 4+ materials, 
this requires Ir-O-Ir bond angles far from 180° (for 180° bonds crucial complex phases in the j — 1/2 wavefunctions 
cancel). Our pyrochlorc model here is an example, as is the two-dimensional honeycomb model considered in Ref. [11]. 
Even when TBI physics is absent, the strong enhancement of spin-orbit physics by correlations must play an important 
role in a great many materials. 
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DEFG02-07ER46452 (D.A.P.). The research facilities at the KITP were supported by the National Science Foundation 
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SPIN-ORBIT EIGENSTATES 



The explicit expressions for the wave functions of spin-orbit cigcnstatcs (j = 1/2 doublet and j — 3/2 quadruplet), 
are given in Eq (1). The spin projections are taken in the axes bound to oxygen octahedra of each Ir site (called local 
octahedral axes from now on). 
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DIMENSIONLESS HOPPING MATRICES 



The hopping matrices are constructed as follows. 

Let us first obtain the hopping matrices for hops between the t 2g states in local octahedral axes. A pair of Ir sites is 
connected via an oxygen site nearest to the two Ir ones. As far as the Ir-0 coupling is concerned, the possible hopping 
matrix elements are restricted by symmetry: if O site lies, say, on the z-axis of the octahedral system of coordinates 
of a particular Ir site, then electrons can only hop from O's p y onto Ir's d yz orbital, or from O's p x onto Ir's d zx 
one. Analogous rules apply to other placements of O sites, and we encode them by introducing the corresponding 
Ir-0 hopping matrices, t^ v z , where the subscript says which axis a particular oxygen resides on, and the superscript 
shows whether it has positive or negative coordinate: Since p-orbitals have odd parity, this changes the sign of the 
hopping matrix element. Explicitly, these matrices, connecting p x ,p y ,p z to d yz ,d zx and d xy orbitals, read 

/ \ / 1 \ / 1 \ 

r± = ±V pd7T 1 ; r± = ±V pd7r ; r f - ±V pd7T 10 0. (2) 
\ 1 / \1 0/ \ 0/ 

Here V pd -n is the common matrix element, parameterizing the strength of the hopping. 

To "integrate out" the oxygens we need to view their p orbitals from the same system of coordinates. We choose 
it to be the global cubic coordinate, and introduce matrices that rotate the cubic axes into the octahedral ones 
of site i: 




The transformation law for p orbitals creation operators is pf n = R z mn Pm , where the index i = ... 3 labels sets 
of the octahedral coordinates, operators with a single index correspond to p orbital in cubic axes, and n,m = 1,2,3 
correspond to x,y,z, respectively. Integrating out the oxygen orbitals, we obtain the hopping matrices between local 
t 2g orbitals: 



t u ' =tr ii ,{R i ) T R i 'r i H, 



(4) 



2 



where t^i is the appropriate choice of r-matrices for hopping from site i onto site i' . For example, for hopping between 
sites and 3 we have T03 = r~ , and T30 = t+. As mentioned in the main text, t oc V 2 d7V /A, where A is the energy 
difference between the Ir d and O p states. 

Ir-0 hopping conserves the electron spin, thus matrices t m are identity in the spin space. However, to project on 
the local spin-orbit eigenstates we would like to view the spins in local coordinate system. To this end we introduce 
matrices which are spinor representations of 0(3) rotations R^: 
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Using these matrices, the spin-dependent dimensionless hopping amplitudes describing hopping of electrons with 
spin projections measured in local axes (index a) between local ti g orbitals (index I) are written as 



rpil' 

± 1„ h 



Tii'(R i ) T R i TiH 



(pW) + D« 



(6) 



The final step consists in changing the basis to that of the spin-orbit eigenstates, listed in Eq. (1). We denote these 



basis states as 



M a j 5 , where index a 



1 ... 6 numbers states from top to bottom in the column in the left 



hand side of Eq. (1), and each set of the coefficients M Qj ; CT is essentially the ath row of the matrix in the right hand 
side of Eq. (1). Using these notations, we finally arrive at the sought dimensionless hopping matrices between the 
spin-orbit eigenstates in local axes: 



11' .era' 



i(R l ) T R l n, 



£)(*') 



(7) 



We will not give the explicit matrices here for space reasons. 



FU AND KANE CRITERION FOR TOPOLOGICAL INSULATORS 



Here we outline the proof that for large X/t the ground state of the system is a strong topological insulator. 
The analysis proceeds in the following way. According to Fu and Kane [1], to extract the topological properties 
of an insulator with inversion symmetry one needs to study the parity eigenvalues of filled bands at time-reversal- 
invariant momenta (TRIM). Below we use the notations of Ref. [1]. Then the eight TRIM are written as Ti =nin2n3 = 
(n\Gi + n-iGi + n^G^)^, where ^2,3 = 0, 1, and Go, 3 arc the reciprocal lattice vectors (i.e. those of the BCC 
lattice, as the direct one is the FCC in our case). 

If the parities, ^ m (Ti), of the occupied Dc mds (index m) are known, then one finds the Z2 topological invariant v 
that distinguishes the strong topological phase from the trivial and weak topological ones as 

2 

= n n ^(ro. (8) 

l m—l 

Here we took into account that for X/t — > 00 there are eight relevant bands, four of which are filled, so the index m 
takes on values 1,2. 

If we choose site as the inversion center (see Fig. 2 of the main text for details), parity acts on a spin-orbit 
cigenstate centered at basis site i of Bravais lattice site R as P4>i{R) = ipi(—R—2bi), where hi are the basis vectors: b Q = 
a(000), 6x = a(011)/4, b 2 = a(101)/4, 63 = a(110)/4, a is the FCC cubic cell lattice parameter. In A:-space at TRIM the 
operator of inversion is a diagonal matrix: P(Ti) = diag(e 2ir ' 6 °, e 2<r,6 °, e 2tTlb \ e 2lTlb \ e 2ir,l>2 , e 2ir ' b2 , e 2lFlb3 , e 2ir ' b3 ). 
Using this fact and diagonalizing the noninteracting part of the hopping Hamiltonian (see Eq.(l) in the main text), 
we can easily find the eigenvalues of P in the subspace of wave functions corresponding to the occupied bands. The 
distribution of ±1 among £ m (r,) depends on the choice of the inversion center, however their product determining 
v, Eq. (8), is not, and we obtain v = 1. Therefore, we conclude that we deal with a strong topological insulator for 
large enough A. 

The other three ("weak") Z2 topological invariants, (^1,^2,^3), can also be obtained knowing £ m 's. While v is 
given by a product of evaluated at all TRIM, each of is determined by a product over those TRIM that reside 



3 



in the same plane [1]: 

2 

(-1)^= LI II &rn{Tl=n m ). (9) 

rii — l,nfc^i— 0,1 m— 1 

An explicit calculation gives v\ = v 2 = ^3 = 0. 



SOLUTION OF THE SELFCONSISTENT MFT 

Following Florens and Georges, we substitute e tf>Ri = xm + iym, with the constraint x Ri + y Ri = 1, again 
treated on average by introducing the corresponding Lagrange multiplier, p. As explained in the main text, we 
set the Lagrange multiplier h = 0, which makes x Ri and y Ri independent. The rotor Hamiltonian, Hg, re- 
duces to a set of coupled Harmonic oscillators on the pyrochlorc lattice. It is easily diagonalizcd: each of x Ri 
and ym leads to four bands of "phonons" that have dispersions Wj(k) = \JU\p + t\Qg\wj(k)], where wi^(k) = 

—2 ^1 ± yjl + J2fj,<v cos(fc^/2) cos(fc„/2)^ , and w^^k) = 2, such that the former two are dispersive, and the latter 

are dispersionless. wi(k) is the lowest band with minwi(fc) = —6. Note also that U had to be rescaled, U — ► U/2, 
in order to keep the correct atomic limit [2]. Having solved the spinon and rotor sectors, we can determine the phase 
boundaries in a self-consistent manner. 

The Mott transition is characterized by the change from gapped, p > 6t\Qg\, to gapless (and Bose-condensed) 
rotors, p = 6t\Qg\. In the uncondensed phase the equations for p and Qf read, respectively 




U 

4(p + t\Q e \w 3 (k))' 



Qf --2^^ Wik ^(p + AQeH k y (10) 

where N sit es is the number of sites in the Bravais lattice. These two equations imply that along the line U C (X), 
separating the gapped and condensed rotors, Qf stays constant and equal to 

J_ 1 

Qj = TT^ 5f±-^«0.16. (11) 



Having obtained the value of the renormalization factor for the spinon hopping, we then solve the spinon sector to 
find the chemical potential p(X), calculate Qg(X), and finally obtain from the first of Eqs. (10) evaluated at p = 6|Q£| 



that U c (Q° g (\)) - AQl{X) [\ E, fe -J^=) « 19-5Q§(A). 



To find the line Xfs{U) separating spinons with a Fermi surface from those in the topological insulating state, or 
its analog on the Mott-insulating part of the phase diagram, we need to generalize equations for p and Qf to account 
for Bose condensation of rotors in the lowest band, W\(k). Introducing the amplitude of Bose condensation, Z, the 
equations for p and Qf read, respectively: 




Q -' z + Q nlv7M- (12) 

To proceed, we first note that on the mean filed level the spinon sector is solved in the same way for condensed and 
uncondensed rotors, thus we do know that the ratio X/(Qft) stays the same along the line separating, say, spinons with 
a Fermi surface from those in the topological insulating state. And, of course, it is equal to the same ratio for U = 
(i.e. it is 2.8). Second, we observe that the value of Qg depends on this very ratio, X/(Qft), only. Therefore, Qg stays 
the same along a line with a given X/(Qft). We then solve the spinon problem at X/(Qft) = 2.8 to obtain \Qg\ w 0.26. 



Then for a given U we find p and Qf. The latter determines the sought line simply as Xps(U) — 2.8tQf(^j.). The 
line corresponding to the gap closure at X/(Qft) w 3.3 is obtained in exactly the same manner. 
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